Analysis of the Holzmann-Chevallier-Krauth theory for the trapped quasi-two-dimensional Bose 

gas 



O 

o 

(N 



R. N. Bisset, D. Baillie and P. B. Blakie 

Jack Dodd Centre for Quantum Technology, Department of Physics, University of Otago, New Zealand 

(Dated: March 7, 2009) 

We provide an in depth analysis of the theory proposed by Holzmann, Chevallier and Krauth (HCK) [Eu- 
rophys. Lett., 82, 30001 (2008)] for predicting the temperature at which the Berezinskii-Kosterlitz-Thouless 
(BKT) transition to a superfluid state occurs in the harmonically trapped quasi-two-dimensional (2D) Bose gas. 
Their theory is based on a meanfield model of the system density and we show that the HCK predictions change 
appreciably when an improved meanfield theory and identification of the transition point is used. In this analysis 
we develop a consistent theory that provides a lower bound for the BKT transition temperature in the trapped 
quasi-2D Bose gas. 
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I. INTRODUCTION 

Thermal fluctuations destroy long-range coherence in the 
two dimensional (2D) Bose gas fll However, a 

Berezinskii-Kosterlitz-Thouless (BKT) transition IH I3l to a 
quasi-coherent superfluid state is predicted to occur, and was 
observed first in liquid helium thin films in 1978 |i5|]. More 
recently evidence for this BKT transition was reported for a 
dilute Bose gas by the ENS and NIST groups ||1 S S H in 
harmonically trapped quasi-2D systems. 

These experiments raise two important issues that need 
to be dealt with: (i) The trapping potentials are not purely 
two-dimensional for the temperature regimes considered, with 
some thermal excitation in the tight direction, (ii) The weak 
harmonic confinement in the 2D plane of the system intro- 
duces finite size effects, and competition between interaction 
and potential energy of the system. These issues have made 
the low temperature phase diagram of this system the subject 
of much debate |10, 11, 12, 13, 14, 15, 16, 17], with reli- 
able predictions only coming recently from classical field and 
quantum Monte Carlo methods Iil8l[l9i l20ll21|] . 

It is of course desirable to have a simple meanfield de- 
scription of the quasi-2D system ll22i l23ll . However, mean- 
field theories are of limited applicability in the critical region, 
where density fluctuations are strong, and it is well known that 
the 2D critical region is large. However, recently Holzmann, 
Chevallier and Krauth (HCK) [23] made a novel proposal to 
use a high temperature Hartree-Fock meanfield theory to ex- 
trapolate into the lower temperature critical regime. They then 
used this theory to estimate the transition temperature (Tbkt) 
as that where the peak phase space density of the system sat- 
isfies the critical value (see Eq. (|9]l) known for the uniform 
pure-2D Bose gas [12 . 13]. 

In this paper our principle focus is to analyze two aspects of 
the HCK theory: (i) The simplification of the interaction term 
to an averaged value that is insensitive of the particular axial 
mode the atoms occupy, (ii) The use of the total areal density 
to identify the BKT transition. In our analysis we use a more 
complete high temperature meanfield theory that avoids the 
interaction simplification used by HCK. We also show that the 
correct generalization of the pure-2D condition for the BKT 



transition to the quasi-2D system involves the areal density of 
the ground axial mode of the system. Through numerical cal- 
culations we show that our improved treatment of these two 
aspects leads to significant differences in our theoretical pre- 
dictions from those of HCK. We also discuss the main limita- 
tion of meanfield theory extrapolation into the critical regime, 
which indicates that our improvements on the HCK theory 
will provide a lower bound on the temperature for Tbkt in 
the trapped quasi-2D system. 



II. FORMALISM 



We begin with the effective low-energy Hamiltonian for 
ultra-cold bosonic atoms 

J I 2m 2 ■> ^ j 



(1) 



where ^'(x) is the quantum Bose field operator that anni- 
hilates a particle at position x, a is the s-wave scattering 
length and m is the atomic mass. The quasi-2D system we 
consider here is realized when the trapping potential is suf- 
ficiently tight in one direction (which we take to be z) that 

Our interest is in the thermal properties of the quasi-2D 
system when there is no condensate present, a regime for 
which Hartree-Fock theory is appropriate. If interactions 
are small compared to hojz then the Hartree-Fock modes for 
the Hamiltonian ([T]i take the separable form i/jka{x, y, z) = 
fka{x,y)^k{z), where the axial modes ^k{z) are bare har- 
monic oscillator states. In the quasi-2D regime the xy-plane, 
for which we introduce the notation r = (.t, y), can be treated 
semiclassically, eliminating the need to diagonalize for the 
modes fka{x,y). However, the axial modes must be treated 
quantum mechanically, and the Hartree-Fock expression for 
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the areal density of the system in the j-th axial mode is 

1 r ,o. 1 
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where the Hartree-Fock energies are 



1 



(2) 
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+ ^{ujIx'^ +uly'^)+jhLO, (3) 
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H is the chemical potential, and 



gki 



dz\Uz)VMz)Y 



(4) 



describes the interactions between atoms in the k and i axial 
modes. Performing the momentum integration in Eq. Q and 
adding up the axial mode densities gives the total areal density 

1 °° 

n{v) = --^Y.^n[l-e^^ - Vj{v)} IkeT)] , (5) 
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where 



= ^{uJxx'^ +t^yy'^) + jf^z + 2^gkjnk{r), (6) 



fc=0 



is the effective potential for atoms in the j-th axial mode, and 
A = h / y/2jrmkBT is the thermal de Broglie wavelength. We 
solve Eqs. (|2]l and (O self-consistently, i.e. by iterating until 
the solutions converge. 



A. Comparison to theory of Holzmann et al. 

A central concern of our work is to compare our meanfield 
theory, as outlined above, to the meanfield theory used by 
HCK ll23ll . The HCK theory is a simplification of our mean- 
field scheme presented above made by taking the interactions 
to be axial mode independent, i.e., changing the meanfield in- 
teraction term to 



(7) 



fc=0 



This approximation has no rigorous justification, but allows a 
closed-form expression for total density. The average interac- 
tion strength used in the HCK model is given by 



an 



m 



dz[p{z)f 



(8) 



where p{z) is the density of a single atom in a harmonic os- 
cillator of frequency uo^ at temperature T. 
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Figure 1 : Comparison of 2D phase space densities for systems at (a) 
T = 221 nK, (b) T = 172 nK, and (c) T = 150 nK. Our meanfield 
model (solid), HCK model (dotted) and ideal gas (dashed). Ground 
axial mode areal densities are shown in grey in (c). Calculation pa- 
rameters are A'' — 10'* **^Rb atoms with LL}x,y = 27r x 59 Hz and 
= 27r X 2530 Hz. We note that the parameters in (a) are the same 
as those in Fig. 2 of 1 23]. 



B. Results: 2D phase space density 

In Fig. [T] we compare the two meanfield models and the 
ideal gas model (i.e. Eqs. (|2]i and (|3]l with a = 0) for 
three different temperatures. The density profiles shown in 
Fig- H^) at the pure-2D condensation temperature |24|, 
i.e. T^Bc = \/&Nh^ujx'^yl'r:kB lHH. In this regime both 
meanfield models make similar predictions, and have lower 
density than the ideal gas at the trap center due to the effects 
of repulsive interactions. 

The results in Fig. [Tib) are at a colder temperature where 
the ideal gas is almost saturated, (i.e. approximately at the 
quasi-2D condensation temperature, T'^'^c < T^^q, where 
the central density begins to diverge). Interaction effects play 
a more significant role here, and prevent the density from 
spiking in both meanfield theories. In this regime, and at 
lower temperatures [Fig. [Uc)], the differences between the 
meanfield theories become clearly apparent. This discrepancy 
arises from how interactions are treated in the theories [see 
Eq. (|7]i] in two ways: 

(i) Averaged interaction parameter: The interaction param- 
eter, (7h, in the HCK theory assumes that the various modes 
follow a Boltzmann distribution. As the system becomes de- 
generate, uqX^ ^ 1, this approximation is inaccurate as it 
fails to account for quantum statistical effects that increase 
the ground mode occupation. A comparison of both meanfield 
theories and the Boltzmann prediction of the ground band oc- 
cupation are given in Table U and reveals the increasing dif- 
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add or subtract meanfield contributions for the high energy 
modes which modifies the total density and the parameter C 
used in Eq. (|9]l. The value C ~ 380 only applies to the con- 
tinuous Bose gas. 

The uniform criterion ^ can be applied to the quasi-2D 
trapped gas using a local density approximation. Within such 
an approximation the transition will be spatially dependent, 
and will occur first at the centre of the trap where the density 
is highest. Some care has to be taken in extending the pure-2D 
theory to the quasi-2D case, since the additional axial modes 
correspond to a different meanfield theory. The correct quasi- 
2D extension of result ^ should be applied to the subsystem 
consisting of atoms occupying the ground axial mode, with 
relevant interaction parameter goo^ i-e. 

«"" = ^lnf^), (10) 

where C « 380. 



1. Meanfield validity near the transition 

We note that when condition ( fTOl i is satisfied the meanfield 
theory must be inapplicable as a complete description. How- 
ever, the meanfield theory itself predicts no transition to BKT 
(or EEC) phases and varies smoothly through this tempera- 
ture range. Thus the basis of the HCK approach, the use of 
meanfield theory to estimate the system phase space density, 
would seem to be a reasonable starting point for estimating 
Tbkt in the trapped system. The main deficiency of the 
meanfield approach is that it assumes Gaussian fluctuations, 
whereas near the transition these are suppressed via the for- 
mation of a quasi-condensate iT2l I21I1 . In such regimes the 
Hartree-Fock treatment overestimates the energetic cost due 
to interactions, and predicts the system to spread out more 
and hence have a lower central density. We hence conclude 
that near the BKT transition our meanfield density will be less 
than the actual system density. As discussed in Sec. IIIBI the 
simplifying assumption of the HCK theory leads to a lower 
self interaction for the ground mode which increases the den- 
sity, though for reasons unrelated to quasi-condensation. Thus 
the HCK theory cannot be assured to provide a lower bound 
for the system density. 



B. Using meanfield tlieory to predict the BKT transition 

While both Hartree-Fock meanfield theories discussed fail 
to predict any transition in the quasi- 2D Bose gas, their pre- 
dictions for the density can be extrapolated to estimate 
when the quasi-2D condition for the BKT transition is ful- 
filled (fTOb . This approach was formulated in li23|] . but using 
the total areal density, n, and the averaged interaction param- 
eter, gn, in expression (|9]l. These choices are in contrast to 
the correct quasi-2D extension of the universal result given in 
Eq. (O. 



Table I: Comparison of parameters and theoretical predictions for 
the quasi-2D system considered in Fig. [T] A'^ — 10* *^Rb atoms. 
Interaction parameters, values of the ground axial mode occupation 
iVo = /d^rno(r) for our theory (TVo), HCK theory {N';^^^^), 
and the Boltzmann case (A'^j)^"'*^'). The peak (areal) total phase 
space and ground axial mode phase space densities of our theory 
are shown. In our theory the first few interaction parameters are: 
{900, 901, gu} = (ftV"i){0-1299, 0.0650, 0.0974}. 

ference between the meanfield and Boltzmann results as the 
phase space density increases. 

(ii) Mode independence of interaction: The marked dif- 
ference between the ideal and meanfield solutions [e.g. see 
Fig. [lib)] arises from interaction effects. These effects are 
dominated by the atoms in the ground axial mode. Because 
.9h < <7oo, our ground mode atoms are more strongly inter- 
acting than those in the HCK theory, thus our theory predicts 
no(r) to be lower and broader [see Fig.fTJc)], with less atoms 
in the ground axial mode [see Table H). 

Finally, we comment on validity aspects of our Hartree- 
Fock approach (also see Sec. IIII A lb . First, the separabil- 
ity of our Hartree-Fock description into planar (/^^(r)) and 
axial (harmonic oscillator, ^a;(^)) modes requires that inter- 
action effects are small compared to the axial energy scale, 
i.e. no(0)goo/^'^2 ^ 1- For the results presented in Fig. [1] 
the value of this ratio varies from 0.034 (Fig. [Ha)) to 0.16 
(Fig.Ec)). 



III. BKT TRANSITION IN THE TRAPPED SYSTEM 

A. Monte Carlo analysis of the BKT transition 

The BKT transition occurs when the superfluid density, 
nsF, satisfies nspA^ — 4 [4]. Monte Carlo calculations by 
Prokof'ev et al., 1 12] have characterized the uniform 2D Bose 
gas, showing that this condition, in terms of the total density, 
is 

""i^-^lnf— ), (9) 

where C = 380 ± 3, g is the 2D interaction strengfli ifHfTsll . 
and the subscript 2D emphasizes that this result is for the 
pure-2D Bose gas. An important point made in Ref. [12] is 
that the long wavelength behavior of all 2D weakly interacting 
IV'I* models is universal at the transition point. However, dif- 
ferences between models emerge in high energy modes, which 
contribute to the total density at the critical point, but do not 
affect the strongly fluctuating critical region. Hence, we can 
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Figure 2: Predictions for thie BKT transition temperature as a func- 
tion of A'^. Our theory (squares), HCK tlieory (triangles), quasi-2D 
BEC (points) and pure 2D BEC temperature (crosses). All other pa- 
rameters as in Fig.[T] 



C. Comparison of meanfield predictions for the BKT 
Transition 

We now explore the difference in the meanfield- 
extrapolated predictions for the BKT transition temperature 
(Tbkt) between (i) our meanfield theory using Eq. (fTOl i and 
(ii) the HCK theory using n"" = A"^ {h^C/mgn). These 
predictions are shown in Fig. |2] for systems of various total 
particle number. Over the entire range of atom numbers con- 
sidered our Tbkt estimates are significantly lower than those 
of the HCK theory, indeed, for the system with 8 x 10'' atoms, 
our prediction is approximately 50 nK lower. However, we 
note that for the largest atom numbers considered, the ratio of 
^bTekt /fi^z is sufficiently large that the system is crossing 
over to the 3D regime. 

A basic comparison of the two Hartree-Fock theories has 
akeady been discussed in Sec. IIIIBI however two main fac- 
tors are responsible for the appreciable differences in their 
extrapolated-predictions of Tbkt'- 

(i) Smaller interaction parameter: The interaction pa- 
rameter used by HCK to determine the critical density (i.e. ^h) 
is smaller than ours (500) [see TableJ). Noting that the interac- 
tion parameter only effects the critical density logarithmically 
([Toll), this difference leads to a negative shift in the HCK pre- 
diction for the critical temperature. 



(ii) Use of meanfield density: Using the total areal density 
to judge when the system is critical causes the HCK theory to 
predict a higher critical temperature than our approach which 
instead uses no(r). Furthermore the HCK meanfield theory 
predicts higher densities through the use of the averaged in- 
teraction parameter, as discussed earher. A comparison of 
the various areal densities of these two theories in the criti- 
cal regime for N — 10* atoms is given in Fig.fTfc). 

Effect (ii) dominates for the results in Fig. |2] increasingly 
so for larger numbers of atoms. 



IV. CONCLUSIONS 

In this paper we have outlined a systematic Hartree-Fock 
meanfield theory for the trapped quasi-2D Bose gas, which 
should provide a good description above the BKT critical re- 
gion. A central concern of this paper has been to compare our 
theory against that of HCK [23], which involves an unjustified 
simplification of the meanfield interaction term. We show that 
the density profiles predicted by these theories disagree as the 
temperature decreases. 

We have also considered how to extend the BKT critical 
density condition to the trapped quasi-2D Bose gas. Using 
this extension we are able to extrapolate our meanfield theory 
to predict the critical temperature. We show that the HCK 
extrapolation, which uses a different local density condition to 
identify the BKT critical point, predicts a significantly higher 
value for Tbkt- 

Finally, we would like to emphasize that our treatment for 
extrapolating meanfield theory to predict the BKT transition 
will likely under-estimate Tbkt- This arises because the as- 
sumed Gaussian fluctuations of the Hartree-Fock theory re- 
duce the system density compared to that of the actual system, 
which forms a quasi-condensate as a precursor to the BKT 
transition. Thus our theory will provide a lower bound on 
where the system will obtain the critical density expected for 
the BKT transition to occur. 
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